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Abstract 

We propose a set of photonic crystals that realize a nonlinear quantum Rabi model equivalent 
to a two-level system driven by the phase of a quantized electromagnetic field. The crystals are 
exactly soluble in the weak-coupling regime and their dispersion relation is discrete. The system is 
diagonalized by normal modes equivalent to a dressed state basis. In the strong-coupling regime, 
we use perturbation theory and find that the dispersion relation is continuous and give the normal 
modes of the crystal in terms of continued fractions. We show that these photonic crystals allow 
state reconstruction in the form of coherent oscillations in the weak-coupling regime. 
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I. INTRODUCTION 



Photonic crystals as classical simulators of quantum processes have been the focus of 
attention in recent years.'^^ In particular, it has been shown that the so-called quantum 
Rabi model describing the interaction of a two-level system with a quantum field may be 
realized by photonic superlattices.'^ The quantum Rabi model in the weak-coupling regime, 
i.e. the Jaynes-Cummings model,'^ describes a variety of quantum mechanical systems 
that have been experimentally implemented; e.g. cavity-quantum electrodynamics (cavity- 
QED)^^, ion trap^ and circuit-QED.'^^ Strong-coupling is not feasible in a majority of 
simple quantum optical systems but photonic crystals provide a classical realization of the 
quantum model in all coupling regimes.'^ 

In quantum optics, diverse non-linear models describing the interaction between a two- 
level system and a quantum field have been proposed as deformations of the Jaynes- 
Cummings model.'^^'^ One example of these nonlinear models is the Buck-Sukumar (BS) 
model where the atom-field coupling depends on the intensity of the quantum field.'^ The 
BS model, which is exactly soluble and does not have a feasible experimental representation, 
unless it is classically realized in a couple of binary photonic crystals where the coupling 
depends linearly on the position of the waveguide, helps in understanding the apparition of 
collapses and revivals of the two-level inversion in the radiation-mater systems. 

In the following, we propose a semi-infinite photonic crystal realization of a non-linear 
model describing a system where the atom-field coupling depends only on the phase of the 
quantum field; something that is missing in the literature up to our knowledge. Then, we 
find the exact dispersion curves and normal modes of the waveguide lattice in the weak- 
coupling regime. In the strong-coupling regime, the dispersion relation is continuous and we 
find the normal modes as continued fractions. The transition from discrete to continuous 
spectrum, appearing in our photonic crystal, does not show in the spectra of RabP^'^ and 
B^^ models which are always discrete. Thus, parameter sets delivering coherent oscillations 
in Rabi or BS models only produce coherent oscillations in the weak-coupling regime of our 
model. 
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II. THE MODEL AND ITS PHOTONIC CRYSTAL ANALOGUE 



Let us consider the Hamiltonian describing a two-level system driven by just the phase 
of a quantum field, 



H = Ufd^a + -jd, + A [e'* + e-^Pj a^, (1) 

where the exponential of the quantum phase operator is given by the Susskind-Glogower 
operatoi'^ 

e'^ = V = -i=a. (2) 

Where the field mode of frequency is described by the annihilation(creation) operators 
a(d^), the two-level system of transition frequency uq by Pauli matrices CTx,y,z, and their 
interaction by the real coupling A. It is possible to separate this system in two uncoupled 
Hamiltonians, 

= ujn T y (-1)" + A (fi + fit) , (3) 
belonging to one of two parity chain basis, 

|+,n) = 5tn|o,^), (4) 
|-,n) = 5t"|0,e), (5) 



defined such that parity. 



n = -a,(-if, (6) 



is conserved, (±,n|n|±,n) = ±; the bases annihilation(creation) operator is given by B 
VaxiB"^ = V^^cT^) and the number operator is defined as n|±,m) = m|±,m). 
By defining the general state. 



oo 

'(±)l 



j=0 

the equations of motion for any given initial state under the dynamics given by Hamiltonian 
([T]) are reduced to the differential set 



id,Sf> 



where the shorthand notation dt has been used for the partial derivative with respect to t. 
This differential set is equivalent, up to a phase and substituting t — )■ z, to that describing 
the propagation equation of a classical field through a photonic waveguide lattice. In this 
equivalent photonic waveguide lattice, £j is the amplitude of the field at the jth waveguide, 
the waveguides are homogeneously coupled, and the refraction indices grow proportional to 
ujf and to their position on the lattice plus a position depending bias proportional to ujq/2. 

For the sake of simplicity, hereby we will refer to the photonic crystals as iJ+ or 
depending on the sign of Eq.(|8|. In order to construct these photonic crystals, one can 
choose to either implement a static. Fig. [T]^a), or dynamic. Fig. [l]^b), relation between 
parameters Uf and wq; bending the waveguides along a circle introduces an index gradient 
inversely proportional to the wavelenght of the impinging light .'^ 




FIG. 1. (Color online) Two schemes to produce the set of two photonic crystals realizing the 
phase driven two-level atom. A semi-infinite set of homogeneously coupled waveguides where the 
refraction index behaves as the function n^-^^ oc ujfj =F ^(—1)-'. (a) Straight waveguides produce a 
static relation between parameters ujf and ojq. (b) Circularly bent waveguides produce a response 
to the incident frequency of light producing a dynamic relation between loj and ujq. 



III. DISPERSION RELATION 

Up to our knowledge and means, it is not possible to find an exact dispersion relation 
for the photonic crystals described above but it is possible to separate the relevant coupling 
parameter in two regimes, weak and strong, in order to obtain some results. In the first of 
these regimes, we can borrow techniques from quantum optics and find an exact dispersion 
relation. While on the last, we can only deal with the problem through perturbation theory. 
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A. Weak-coupling regime: \<^uj,ujo. 



In the weak coupling regime one can find the exact spectrum and normal modes of each 
one of the photonic crystals described by the differential set ([s]) by taking a step back and 
implementing the rotating wave approximation in the Hamiltonian of the system, 



HnwA = ujfa^a + ^a, + A (^e^-^a^ + e~'^a_j , (9) 

before establishing the classical analogue. Then, it is simpler to find the spectrum and 
normal modes in this representation by using the basis set {\n,e), \n + l,g)} belonging to 
the manifold with n + 1 excitations. This leads to the discrete spectrum, 

i?±,n = cu/(ri + ^^ n = V5' + A\^ (10) 

where the proper states are given by 

\n,±) = a±\n,e) + P±\n + 1, g), (11) 



with 

a± -5±n 



(12) 



(3± 2A 

Note that \0,g) is an eigenstate of the Hamiltonian with energy -E-,o = — i^o/2. 

In our photonic crystals, Eq. (jsj), we can approximate the dispersion relation, equivalent 
to the discrete spectrum found above, by proposing a collective proper mode and realizing 
that the three-term recurrence relations can be summarized by the tridiagonal matrix, 

H±,w = Hlt^ + P, (13) 

'^fjT^{-iy]kj, (14) 
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in, = ^ + '^^+1 J , (15) 

where the notation {M)-j stands for the (i,j)th term of Matrix M and the symbol 6a,b is 
Kronecker's delta. As A ^ ujf,u)o, we can treat matrix P as a perturbation on matrix Ho 
and find the eigenvalues of H± up to second order corrections as the first order correction is 
equal to zero. Thus, we obtain the approximated dispersion relation, 

uigY^^^u,gT^{-iy(l + -^]. (16) 
Figure [2] shows good agreement between the dispersion relation given by the exact eigenvalues 



in the rotating wave approximation, Eq. (10), and the perturbation approach, Eq. (16), at 
zero and second order. 
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FIG. 2. (Color online) A segment of the dispersion relation in the weak-coupling regime. Exact 
closed form from rotating wave approximation (solid black), 0th order perturbation (dashed blue) 
and second order perturbation (dotted red) are shown. We have used a detunning given by ujq = 
1.1 u;f. 

B. Strong-coupling regime: A^WjWq- 

In the case when the coupling parameter is larger than the field and transition frequencies, 
sometimes also called deep-coupling regime, it is possible to write the three term recurrence 

as 



^0/ , 



Si J. 



(17) 
(18) 
(19) 



Notice that lead order matrix is diagonalized via Hq = VAV~^, where the diagonal matrix 
A contains the values of the dispersion relation in its diagonal and the matrix V has as 
columns the coefficients of the normal modes given by 



•1 2\ J 



(20) 



where the function ?7„(x) is the nth Chebyshev polynomial of the second kind evaluated at 
X] i.e., the dispersion relation for this case is continuous. 

At this point it is irrelevant, but first order correction for the dispersion relation delivers. 



uj{q) 



± 



/x(g) + / dfi{q)- 



k=0 



2A 



E 



j=0 



2A 



(21) 
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IV. COLLECTIVE MODES 



For any given set of parameters, decomposition in normal modes delivers the three term 
recurrence mentioned before, 



uj{q) 



-.(±■9) 



uj{q) 



(22) 
(23) 



with 



(±) • _L 



(24) 



where the coefficient cj^''^^ is the A;th coefficient of the qth collective mode corresponding to 



the proper value uj{q). These coefficients are given by. 



j-l(±,q)(±,q) 



-0 



(25) 



where we have used the continued fraction 



A 



uj{q) - 4S - A.gf 
A 



(±,9) ' 



uj{q) - 4+i - 



A2 



(26) 
(27) 
(28) 



Lo(q) — a 



{±) 

J+3 



where one can always set cq — 1 and normalize the semi-infinite set later. Notice that 
in the weak-coupling case a normal mode equivalent to that found in the rotating wave 
approximation treatment is recovered as the continued fraction is cut at the second term 
due to the parameter A^ being negligible compared to the frequency of the field. 

In any given regime, wc can take the continued fraction result and use it to write the 
eigenvectors of the quantum Hamiltonian in the reduced form. 



\e±,r. 



j=0 



(29) 
(30) 



where Cq is chosen such that {e±^m\e±^n) = S„ 



V. PROPAGATION EXAMPLES 



The dressed state basis that diagonahze our model, Eq.(ll), imphes that starting in a 
state of the kind |j, e) or \j + l,g) with j > will produce coherent oscillations. Such 
an initial state translates to laser light impinging the jth or (j + l)th waveguide of one 
of the lattices; the case j even (odd) corresponds to the crystal ). Figure |3]^a) 

shows the propagation of light impinging at the 0th waveguide of our photononic crystal 
H_ in the weak-coupling regime which is equivalent to an initial state |0, e) in the quantum 
optics model. Accordingly, we observe the intensity oscillate between the 0th and first wave 
corresponding to an oscillation between the |0,e) and states in the quantum optics 



model. Figure |3j(c) presents a detail of the intensity at the 0th and first waveguide. This 
normalized intensity at the 0th waveguide is proportional to the probability of finding the 
time evolution of the quantum system back in the initial state, Jq ^ -P-.o, with 

P_,o = |(-,0|V^(t))P, |^(0)) = |-,0). (31) 

In quantum optics literature, it is known that the quantum Rabi model produces coherent 
oscillations in the two-level system inversion for uq = 0.^^ In our model, due to the continuous 
spectra, all proper states are scattering states and at most we observe partial recovery of 
the original state when the field starts localized at a given waveguide. Figure ^h) shows 
the propagation of light impinging at the 0th waveguide of the photonic lattice if+ in 
the strong-coupling regime, which is equivalent to an initial state \0,g) in the quantum 
optics model; Figure focus on the intensity at the first two waveguides. Again, the 
normalized intensity at the 0th waveguide is proportional to the probability of finding the 
time evolution of the quantum system back in the initial state, Jq ^ P+,o = l(+;0|'?/'(t))p, 
with \ip{0)) = |+,0), for this case. 

In both weak- and strong-coupling simulations a finite photonic lattice of size 5000 were 
used to produce the numerical results. 



VI. CONCLUSION 



We have proposed a set of two photonic crystals that classically simulates a new radiation- 
matter interaction where a two-level system is driven by the phase of a quantum field. 
We show that it is possible to determine exactly the dispersion relation of the photonic 
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{wo,A} = {l,0.1}cj/ {tjo.A} = {0,2}cj/ 




FIG. 3. (Color online) Examples of propagation in our positive parity photonic crystal. The left 
column (a,c) shows the case of weak-coupling, A = 0.1 cof, and the right (b,d) the case of strong- 
coupling, X = 2 ojf. The first row (a,b) depicts intensity propagation on the first ten waveguides of 
a total of five thousand when light impinges the zeroth waveguide. The second row (c,d) shows the 
normalized intensity at the 0th (solid black) and first (dashed blue) waveguides. The dimensionless 
time parameter LVjt is equivalent to the typical dimensionless propagation parameter. 

waveguide lattices in the so-called weak-coupling regime and that in the strong-coupling 
regime we can use perturbation theory to approximate the dispersion relation up to second 
order perturbation. In the first case, the dispersion relation is discrete and, in the latter, 
continuous. The normal modes of the crystals are easily expressed in terms of continued 
fractions as a function of the dispersion relation. Our optical realization of a phase driven 
two-level system provides a scheme to explore an interesting process that is not accessible 
by usual means; i.e. cavity-QED or trapped ions. 
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